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A  collisional  radiative  model  is  presented  for  analyzing  the  xenon-propelled  Hall  thruster  optical 
radiation  based  on  apparent  electron  and  ion-impact  emission  cross  sections  associated  with  lines  in 
the  visible  and  near-infrared  region  of  the  spectrum.  The  emission  cross  sections  of  selected 
near-infrared  emission  lines  are  incorporated  in  a  collisional-radiative  model.  The  effect  of  stepwise 
excitation  via  metastable  states  on  the  derived  line  intensities  for  emissions  from  Xel  5p56p{6p’) 
levels  is  evaluated.  Meanwhile,  visible  Xell  emissions  are  shown  to  provide  plasma  densities  at  high 
electron  temperature  conditions.  The  electron  temperature  and  spatial  ion  number  density 
distribution  were  determined  from  the  luminescence  spectra  measured  in  the  discharge  and  plume 
near-field  plasma  of  the  Hall  thruster,  the  TSNIIMASH  D-55  anode  layer  thruster.  The  results  are  in 
good  agreement  with  the  probe  measurement  data  and  simulations  reported  in  the  literature  for  the 
same  thruster.  The  analysis  of  the  Hall  thruster  Xel  near-infrared  spectra  demonstrates  that  the 
neglect  of  ion-atom  collisions  results  in  an  erroneous  electron  temperature  determination  at  electron 
temperatures  below  10  eV.  ©  2006  American  Institute  of  Physics.  [DOI:  10.1063/1.2195019] 


I.  INTRODUCTION 

The  extraction  of  plasma  microparameters  from  spectral 
line  intensities  observed  in  the  radiation  emanating  from  xe¬ 
non  propelled  Hall  effect  thruster  (HET)  plasmas  requires  a 
collisional  radiative  model  (CRM)  that  takes  all  relevant  ex¬ 
citation  and  deexcitation  mechanisms  into  account.  While  in 
many  cases  the  assumption  of  equilibrated  electron  energy 
distributions  is  justifiable,  this  cannot  be  stated  for  the  atom 
and  ion  electronic  state  distributions,  as  previously 
reported.1-3  For  strong  emission  lines  in  the  rarefied  condi¬ 
tions  of  the  HET  plasma,  the  corona  radiation  model4  may  be 
considered  as  a  starting  point  for  an  approximate  CRM.  The 
model  assumes  population  of  a  particular  state  through  elec¬ 
tronic  collisions  and  cascade  radiative  transitions  from  upper 
excited  states,  while  depopulation  of  all  excited  states  goes 
through  radiative  decay.  In  many  earlier  HET  optical  diag¬ 
nostic  studies  a  corona  model  was  applied  or  assumed.  Only 
few  efforts  are  known  to  us  that  implemented  a  more  com¬ 
plex  radiation  model  to  analyze  a  HET  xenon  plasma.5,6 
Among  the  corona  model  studies,  Meezan  et  al?  are  the  sole 
authors  that  discussed  the  appropriateness  of  the  model  at  the 
particular  conditions. 

Due  to  the  complexity  of  the  processes  involved  in  HET 
xenon  plasmas,  the  corona  approximation  usually  falls  short 
in  rendering  many  of  the  plasma  emission  lines.  One  source 
of  the  problem  is  the  very  small  energy  intervals  of  Xel 
excited  levels  such  as  the  np5(n  +  \)p  and  np5nd  states.  Con¬ 
sequently,  low-energy  electrons  can  have  very  high  prob¬ 
abilities  in  inducing  transitions  between  these  states.  Another 
source  of  problems  in  the  corona  model  is  the  possibility  of 
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stepwise  excitation  via  metastable  states,  the  importance  of 
which  has  been  previously  demonstrated  for  xenon.  Bu¬ 
grova  et  al?  also  found  that  the  corona  model  failed  to  de¬ 
scribe  the  intensity  of  some  plasma  emissions  in  the  ioniza¬ 
tion  region  of  a  stationary  plasma  thruster  (SPT).  They 
introduced  a  more  complex  CRM  including  stepwise  excita¬ 
tion,  and  obtained  a  better  agreement  with  experiment. 
Meanwhile,  a  final  important  shortcoming  of  the  corona 
model  is  the  neglect  of  heavy  particle  collisions,  which  have 
recently  been  proposed  to  play  an  important  role  in  the  emis¬ 
sion  excitation  process  in  HET  plasmas  given  the  significant 
cross  sections  at  the  respective  high  ion  energies.9 

A  major  deficiency  of  all  models  has  been  a  lack  of 
experimentally  determined  emission  excitation  cross  sections 
for  all  relevant  Xe  plasma  processes.  In  a  companion  paper 
(Part  I),10  we  report  apparent  emission  cross  sections  associ¬ 
ated  with  electron-  and  ion-Xe  atom  collisions  for  a  set  of 
near-infrared  (NIR)  and  visible  (VIS)  lines.  Since  the  re¬ 
ported  cross  sections  include  cascades  from  states  above  the 
upper  radiating  state,  they  incorporate  the  information  of  a 
significant  fraction  of  the  level-to-level  radiation  kinetics  and 
the  related  level  populations.  In  the  present  paper,  we  dem¬ 
onstrate  how  the  use  of  these  cross  sections  considerably 
simplifies  a  CRM  involving  the  Xel  NIR  lines.  The  proposed 
CRM  includes  a  stepwise  excitation  mechanism  without  an 
explicit  knowledge  of  the  cross  sections  for  metastable  exci¬ 
tation  and  excitation  out  of  metastable  states.  The  proposed 
model  also  clearly  confirms  earlier  assertions9  that  the  inclu¬ 
sion  of  ion  processes  is  relevant  for  properly  determining  an 
electron  temperature  in  plasma  regions  with  relatively  low 
characteristic  electron  temperatures. 

After  discussing  the  salient  features  of  the  model,  it  is 
applied  to  determining  the  electron  temperature  and  charge 
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density  distribution  in  a  TSNIIMASH  D-55  anode  layer 
(TAL)  thruster.11-14  We  conclude  the  paper  with  a  discussion 
on  the  benefits  and  limitations  of  the  present  diagnostic 
method,  and  provide  suggestions  for  further  improvements 
of  the  model. 


II.  COLLISIONAL  RADIATIVE  MODEL 
A.  General  model 

In  this  section,  a  CRM  is  developed  to  permit  quantita¬ 
tive  analysis  of  a  HET  plasma  based  on  passive  optical  mea¬ 
surements.  A  CRM  establishes  the  relationship  between  the 
HET  microparameters,  such  as  particle  number  densities  and 
electron  energies,  and  the  optical  radiation  observed  in  the 
HET  plasma.  As  delineated  in  Part  I,  the  optical  radiation  is 
caused  primarily  by  the  following  collision  processes: 


e  +  Xe(p)  — >  Xe *  +  e  , 

(1) 

— >  Xe',+ *  +  (q+  l)e'\ 

(2) 

e~  +  Xe+(p)  — *  Xe‘,+’  +  qe~. 

(3) 

Xe+  +  Xe  — *  Xe+*  +  Xe*, 

(4) 

Xe2++  Xe  — *  Xe2+*  +  Xe*, 

(5) 

— Xe+*  +  Xe+\ 

(6) 

where  the  asterisks  in  processes  ( 1 )— (6)  signify  excited  spe¬ 
cies  that  either  lead  to  optical  emissions  or  long-lived  meta¬ 
stables  that  transfer  the  excitation  in  subsequent  electronic  or 
atomic  collisions,  q  identifies  the  charge  state,  and  p  identi¬ 
fies  all  long-lived  states,  i.e.,  ground  and  metastable  states. 

The  plasma  radiation  intensity  per  unit  volume  for  a  spe¬ 
cific  emission  line  at  wavelength,  X,  may  be  approximated  by 


A  =  2  ^  [  fAEe)(/  (E,)ucdE, 

4ij-\  [  p  NoJo 

q  "  e  J  0 


(7) 


Ne,  Nq,  Np ,  and  N0  are  number  densities  of  the  electrons, 
ions,  and  xenon  atoms  in  state  p  and  all  xenon  atoms,  respec¬ 
tively.  fe(Ee)  and  fq(Eq)  are  the  electron  and  ion  energy  dis¬ 
tribution  functions,  respectively,  normalized  to  unity.  <rep(Ee) 
is  the  cross  section  for  an  electron  with  energy,  Ee ,  to  excite 
a  Xe  atom  or  ion  from  the  lower  state  p  to  produce  radiation 
at  X.  o$(Eq)  is  the  cross  section  for  an  ion  with  charge  q  and 
energy,  Er  to  excite  an  atom  to  produce  radiation  at  X.  ue 
and  uq  are  the  electron  and  ion  velocities,  respectively. 

Equation  (7)  assumes  that  the  states,  i,  excited  by  elec¬ 
tron  and  ion  collisions,  are  depopulated  through  radiation 
decay  only.  Clearly,  this  assumption  is  not  justified  for  some 
xenon  atomic  and  ionic  states.  The  radiation  decay  condition 
may  be  expressed  algebraically  as 


Trad,i*'dex,i <  1  or  A;  =  2  An t  >  vdcx  i,  k  <  i,  (8) 

k 

where  rrad ,  is  the  radiative  lifetime  of  the  upper  zth  level, 
^dex,/  *s  collisional  quenching  frequency  of  the  excited  zth 
level,  and  Aik  is  the  Einstein  coefficient  for  the  /— >k  sponta¬ 
neous  emission  transition.  Although  all  plasma  particles  can 
quench  a  state,  electron  collision  frequencies  are  signifi¬ 
cantly  higher  than  heavy  particle  collision  frequencies  at 
HET  plasma  conditions.  Consequently,  condition  (8)  may  be 
rewritten  as 


A,-  >  Ne  fe(Eg)ade(Ee)ucdEe  =  Nekde ,  (9) 

J  o 

where  (rdg(Ee)  is  the  electron  collision  depopulation  cross 
section  for  the  zth  level,  and  kde  is  the  associated  rate  coef¬ 
ficient.  The  electron  collision  depopulation  rate  coefficients 
in  xenon  may  reach  very  high  values  characteristic  for 
collisional  transitions  between  states  coupled  by  a  dipole- 
allowed  transition.  Estimates  of  the  cross  sections  and  rate 
coefficients  based  on  the  generalized  Bom  approximation15 
for  the  most  intense  NIR  5pf'6pj-+5p55dj,lsj  transitions  re¬ 
sult  in  cross  sections  of  the  order  of  10_l5-10~14  enr  at  Te 
^2  eV.16  At  an  electron  temperature  of  Te^5  eV,  and  a 
plasma  number  density  Ne^  1011  cm-3,  which  is  typical  for 
the  medium  power  (0.5-5  kW)  HET  plume,1417  condition 
(9)  becomes 

Aj>  2X104s-!.  (10) 

In  the  very  near  field  of  the  plume  and  in  the  discharge 
region,  the  electron  number  density  can  exceed  1012  cm"3 
and  a  large  fraction  of  the  electrons  has  energies  Ec 
^7-10  eV.  At  these  conditions,  significantly  higher  sponta¬ 
neous  emission  rates  are  needed  to  satisfy  the  assumption  of 
pure  radiative  decay, 

A,  >  3  X  106  s-1 .  (11) 

Even  though  the  reported  transition  probabilities  and  excited 
state  lifetimes  in  xenon  differ  significantly,18  23  it  can  be  con¬ 
cluded  from  examination  of  the  upper  levels  of  all  prominent 
Xel  and  Xell  emission  lines  in  the  HET  plasma  that  they 
meet  the  condition  (10)  and,  therefore,  Eq.  (7)  applies  to 
these  lines  in  the  HET  plume.  Since  the  strong  Xell  emission 
lines  have  very  short  radiative  lifetimes,  Eq.  (7)  can  be  ap¬ 
plied  to  analyze  the  ionic  emissions  both  in  the  HET  plume 
and  discharge.  The  list  of  Xel  lines  that  could  be  described 
by  Eq.  (7)  in  the  HET  discharge  and  very  near-field  plume 
plasma  is  fairly  limited.  The  prominent  NIR  lines  associated 
with  5p56p(p')  upper  levels,  for  which  emission  cross  sec¬ 
tions  are  reported  in  Part  I  of  this  article,  are  among  the  few 
Xel  lines  that  meet  condition  (11). 

The  integrals  in  Eq.  (7)  represent  the  emission  excitation 
rate  coefficients.  Equation  (7)  can  be  rewritten  in  terms  of 
the  excitation  rates, 
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where  k^p  is  the  emission  excitation  rate  coefficient  for  elec¬ 
tron  collisions  with  particle  in  state  p\  k\  and  k\  are  the 
emission  excitation  rate  coefficients  for  collisions  of  Xe+  and 
Xe2+  with  neutral  xenon  atoms,  respectively,  and  Nr  Nlt  and 
N2  represent  corresponding  number  densities.  Given  the  low 
densities  of  higher  charged  ions  in  HET  plasmas,  only  the 
singly  and  doubly  charged  xenon  ions  are  taken  into 
consideration. “4~26  Equation  (12)  neglects  the  excited  ionic 
states  and  collisions  of  ions  with  metastable  atoms.  Presently, 
it  is  assumed  that  the  respective  rate  coefficients  are  compa¬ 
rable  with  the  respective  ground  state  rate  coefficients. 

The  excitation  rate  coefficients  may  be  expressed  di¬ 
rectly  through  corresponding  cross  sections  if  the  electron 
and  ion  energy  distribution  functions  are  specified.  In  the 
HET  plasma,  a  Maxwellian  distribution  function  is  most 
commonly  assumed  for  electrons.  Ion  energy  distributions  in 
the  HET  plume  are  featured  by  several  peaks  of  different 
heights  and  widths.24-26  Each  peak  represents  a  group  of  ions 
of  different  charge.  The  singly  charged  ions  are  the  most 
abundant  with  a  peak  typically  centered  at  an  energy  some 
tens  of  eV  below  qeV,  where  e  is  the  electron  charge  and  V 
is  the  thruster  anode  voltage.  Doubly  charged  ion  signals  are 
noticeably  weaker  and  located  at  about  twice  the  energy  of 
the  singly  charged  ion  peak.  Widths  of  the  peaks  are  about 
several  tens  of  volts.  In  Part  I  we  report  that  ion  collision 
excitation  cross  sections,  cr^  and  cr^,  vary  smoothly  in  the 
E/q  range  of  100-400  eV,  with  a  slight  incline  with  energy. 
Therefore,  in  a  first  approximation  it  can  be  assumed  that  the 
ionic  cross  sections  are  independent  of  energy  within  the 
range  of  the  ion  beam  energy  distribution  in  the  HET  plume. 
From  these  considerations,  excitation  rate  coefficients  can  be 
determined  from  the  electron  temperature  and  the  ion  beam 
energy  using 


where  me  is  the  rest  mass  of  the  electron  and  M  is  the  atomic 
mass  of  xenon. 

The  number  of  parameters  in  equations  (7)  and  (12)  may 
be  further  reduced  if  plasma  quasineutrality  is  taken  into  ac¬ 
count,  Ne=N[  +  2N2 .  Equation  (7)  then  takes  the  form, 


where  a=N{/Ne  is  the  ratio  of  the  first  ion  (Xe+)  number 
density  to  the  electron  number  density.  Equation  (16)  can  be 
defined  more  narrowly  for  specific  species  emissions.  In  the 


case  of  Xel  emission  lines,  all  Xell  contributions  can  be 
eliminated  from  the  equation,  leading  to 

•/x(Xel)  =  ^Votfe(*So  +  jjrkxem  +«•**  +  , 

(17) 

where  k1^  is  the  rate  coefficient  for  process  (1)  involving 
ground  state  atoms;  k)m  is  the  corresponding  rate  coefficient 
averaged  for  all  metastable  Xe  atoms,  and  Nm  is  the  Xe  atom 
metastable  number  density. 

The  charge-state  parameter,  a,  typically  varies  within  a 
0.7-0.95  range,  depending  on  the  thruster  type  and  operating 
mode. 1,24,25  Relative  abundances  of  the  metastable  atoms, 
Nm/N0,  and  excitation  rate  coefficients,  have  not  been 
determined  at  HET  plasma  conditions.  In  order  to  circum¬ 
vent  this  unknown  parameter,  the  analysis  can  be  limited  to 
lines  that  are  not  coupled  to  metastable  states  through  dipole- 
allowed  transitions.  For  the  NIR  lines  associated  with  the  2p} 
[Paschen  notation  for  5ps6p(p')  set  of  states]  states,  this 
would  include  the  2 p,,  2 p3,  and  2 p5  emission  lines  at  788.7, 
834.7,  and  828.0  nm,  respectively.  These  lines  are  associated 
with  radiative  transitions  to  the  Xel  5ps6s’(\ /2)|  and 
5p56s(3/2)i  (1^2  and  1^4  in  the  Paschen  notation)  states  that 
are  strongly  optically  coupled  to  the  Xe  ground  state  through 
a  vacuum-ultraviolet  transition.  Rate  coefficients  for  excita¬ 
tion  of  2 pj,  2 p3,  and  2 p5  levels  out  of  metastable  states  are 
considerably  smaller  than  those  for  dipole-allowed  transi¬ 
tions,  and  can  be  considered  comparable  with  the  rate  coef¬ 
ficients  out  of  the  ground  state.  This  is  confirmed  by 
recent  work  by  Jung  et  al.  in  studies  of  corresponding  Kr 
levels.27  As  will  become  apparent  from  Eq.  (26)  below,  if  the 
metastable  atom  content  in  the  HET  plasma  is  below  1%,  the 
second  term  in  Eq.  (17)  may  be  ruled  out  when  intensities  of 
Xel  788.7,  834.7,  and  828.0  nm  are  modeled.  A  Xel  CRM 
based  on  these  three-lines  could  thus  be  represented  by 

A(Xel)  =  (<4  +  « '  *i  +  ^*2) •  OS) 

where  all  the  terms  can  be  computed  from  expressions 
(13)— ( 15)  and  the  apparent  emission  cross  sections  reported 
in  Part  I.  The  disadvantage  of  this  model  is  that  it  relies  on 
relatively  weak  lines  at  788.7  and  834.7  nm,  making  reliable 
intensity  measurements  more  difficult.  We,  thus,  consider  it 
useful  to  derive  a  more  general  expression  that  takes  meta¬ 
stable  populations  into  account,  thereby  extending  the  num¬ 
ber  of  lines  available  for  analysis. 


B.  NIR  CRM  including  emissions  coupled 
to  metastable  levels 

The  strongest  NIR  lines  are  associated  with  2 pt  transi¬ 
tions  to  the  5p56^(3/2)2  (IS5)  metastable  state.  Given  the 
high  oscillator  strengths  and  low  excitation  energies  out  of 
this  energy  level  to  the  radiating  5ps6p(2pi)  levels,  rate  co¬ 
efficients,  k)m,  may  be  orders  of  magnitude  higher  than  k 
This  is  fully  confirmed  by  recent  electron-impact  studies  of 
analogous  levels  in  Kr.27  Thus,  excitation  rates  from  states 
coupled  to  this  metastable  state  through  a  dipole-allowed 
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transition  can  be  comparable  to  excitation  rates  from  the 
ground  state,  and  must  be  considered  if  the  more  intense  NIR 
lines  are  to  be  included  in  the  CRM. 

The  steady-state  equilibrium  rate  equation  for  a  particu¬ 
lar  metastable  state  takes  the  form, 

N0(Nekem  +  2  =  Nm(vdm  +  vcm) ,  (19) 

where  kem  is  the  rate  coefficient  for  excitation  of  the  meta¬ 
stable  level  by  electron  impact,  kqm  is  the  rate  coefficient  for 
excitation  of  the  metastable  level  in  a  collision  with  the  first 
(^=1)  and  second  (4=2)  ions,  vdm  is  the  frequency  of  meta¬ 
stable  atom  loss  due  to  diffusion,  and  vcm  is  the  frequency  of 
the  metastable  atom  loss  due  to  collisions.  Mityureva  and 
Smirnov  have  shown  that  in  case  of  excitation  induced  by 
electron-atom  collisions  metastable  states  such  as  the  IJ5 
state  are  predominantly  populated  through  cascade  tran¬ 
sitions.28  Since  the  collision-induced  cascade  rates  are  re¬ 
flected  by  the  apparent  NIR  emission  cross  sections  reported 
in  Part  I,  the  source  term  in  Eq.  (19)  can  be  approximated  by 

N0(Nekem  +  2  =  ¥.2  (*io  +  « •  k\  +  ^*2)  • 

(20) 


whereby  the  sum  on  the  right  side  of  Eq.  (20)  occurs  over  all 
emission  lines  for  which  the  lower  state  is  the  1^5  (J= 2) 
level.  The  selection  rules  for  dipole  transitions  dictate  that 
the  upper  states  of  the  populating  transitions  are  1^7^  3 
states.  Thus,  of  the  eight  NIR  lines,  only  those  associated 
with  upper  states  2/?6(7=2),  2/77(7=l),  2/?8(7=3),  2/79(7=2), 
and  2p[0(J=\)  need  to  be  considered.  We  also  neglect  the 
emissions  from  np  (n> 2)  states,  because  their  intensities 
observed  in  HET  and  beam-induced  spectra  are  much  weaker 
in  comparison  with  those  from  2 pj  states. 

The  wavelengths  of  the  lines  summed  on  the  right  side 
of  Eq.  (20)  are  X  =  823.3,  840.9,  881.9,  904.5,  and  980.0  nm. 
Note  that  the  840.9  nm  line  involving  the  2pn  level  is  asso¬ 
ciated  with  a  very  small  branching  ratio  and,  due  to  its  weak 
intensity,  is  not  recommended  for  diagnostic  purposes.  The 
respective  840.9  nm  emission  cross  sections  can  be  obtained 
from  the  916.3  nm  2p1  line  cross  sections  by  applying  the 
reported  intensity  ratio,  Ig\Jh4\  of  9.9. 29,30  The  lower  state 
of  the  916.3  nm  transition  is  the  \s4  state  that  is  optically 
coupled  to  the  ground  state. 

The  accuracy  of  the  total  rate  coefficient  evaluation  can 
be  improved  by  representing  the  excitation  rate  coefficients 
in  terms  of  relative  units  of  one  line  in  relation  to  all  others 
coupled  to  the  subject  metastable  state.  For  the  evaluation  of 
a  Xel  2 p  emission  line  intensity  of  wavelength,  X,  we  intro¬ 
duce  the  function  Kk(Te,a ),  which  represents  the  ratio  of  the 
total  excitation  rate  of  the  \s5  state  through  emissions  from 
the  2 p  states  to  the  excitation  rate  through  X, 


K\Te,a) 


+  a.*‘+[(l-a)/2]*'} 

i 

*Jb  + a- *}  +  [(! -«)/2]*2 


(21) 


The  sum  in  the  numerator  is  the  same  as  in  Eq.  (20). 
The  relevant  levels  for  the  \s5  coupled  system  are  shown 


2Pm< 


Xel 


FIG.  1.  Simplified  level  diagram  of  2pt  state  system  based  on  the  2 ex¬ 
ample.  lp0:  ground  state;  U5:  metastable  state;  and  1  sA\  optically  coupled  to 
ground  state.  The  figure  identifies  the  transitions  that  play  an  important  role 
in  the  stepwise  excitation  model. 


representatively  in  Fig.  1.  Use  of  Eq.  (21)  then  transforms 
Eq.  .(20)  to 

=  WV.(*Jb  +“•*!  +  ^*2)  Kk(Te . «)  •  (22) 

In  the  loss  term  in  Eq.  (19),  it  can  be  safely  concluded  that 
collisional  depopulation  of  the  metastable  states  is  dominated 
by  electron  collisions.  From  the  known  Xe  atom  quenching 
rate  coefficients31  and  due  to  the  relatively  small  atom  ve¬ 
locities,  it  can  be  concluded  that  the  Xe  atom  quenching  rate 
is  much  smaller  than  the  electron  collision  quenching  rate.  At 
sufficiently  high  electron  temperatures  ( Te>  1  eV),  electron- 
induced  transitions  to  an  upper  2p ,•  level  are  expected  to  be 
significantly  more  likely  than  superelastic  collisions.32  Be¬ 
cause  the  rate  coefficients  for  transitions  between  optically 
coupled  levels  are  much  higher  than  those  for  uncoupled 
levels,  the  collisional  deexcitation  rate  ucm  in  Eq.  (19)  may 
be  given  by  a  sum  of  the  excitation  rate  coefficients,  Th 
corresponding  to  five  \s5— >2p£J=  1 , 1 ,2,2,3]  transitions, 

vcm  =  N^PrTh  (23) 

I 

where  Pr  are  the  branching  coefficients  for  a  radiative  tran¬ 
sition  from  the  upper  2 p,  state  into  the  Xel  ls4  state  that  is 
optically  coupled  to  the  ground  state  (see  Fig.  1).  We  assume 
that  an  atom  excited  to  a  2/7,  level  has  only  two  radiative 
transition  pathways,  either  to  the  Lv4  state  with  probability 
Pr  or  to  the  \s5  state  with  probability  Pm .  All  other  depopu¬ 
lation  processes  are  neglected,  so  that  Pr  +  Pm  =  1  for  all  2/7, 
levels. 

Diffusion  losses  can  be  estimated  from  the  thermal  ve¬ 
locity,  unv  of  the  metastable  atom  and  the  scale  length  of  the 
investigated  plasma,  L, 


vdm  585  “JL-  (24) 

From  the  estimated  excitation  rate  coefficients,  kkm,  based 
on  the  Bom  approximation15  and  the  recent  work  by  Jung 
et  al .,27  and  estimated  diffusion  loss  frequencies,  based  on 
the  direct  simulation  Monte  Carlo/particle-in-cell  (DSMC/ 
PIC)  calculations,17  it  can  be  concluded  that  the  diffusion 
term  can  be  neglected  at  typical  HET  conditions  Ne 
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>  1011  cm  3,  Te> 3  eV  in  the  discharge  and  Ne>  1010  cm  3, 
Te>  1.8  eV  in  the  plume, 

* ’em  =  Nr^  PrTi  >  •'dm ■  (25) 

i 


Insertion  of  Eqs.  (22)  and  (23)  in  Eq.  (19)  results  in  an  ex¬ 
pression  for  the  relative  metastable  content,  Nm/N0, 


Nm  _{!&+ a  ■$  +  [(  1  -  a)/2]^}Arx(r<t  a) 

"o  2  PrTi 

i 


(26) 


At  these  conditions,  the  metastable  content  is  determined  to 
be  in  the  range  of  0.01%-0.3%,  depending  on  Te,  a  and  the 
thruster  body  scale. 

Insertion  of  Eq.  (26)  in  Eq.  (17)  provides  the  following 
expression  for  the  emission  intensity  for  a  particular  line  of 
wavelength  X: 


he 


MXtl)  =  —N0Ne[kxe0+ak\  + 


1  —  a 


47 t\ 

X 

=  ^N0Ne{k^+ak\  + 


S  PrTi 


1  -  a 


1  + 


K\Te,a) 


(27) 


The  term  Kx(Te,a)/ £x  in  Eq.  (27)  represents  the  fractional 
increase  in  intensity  due  to  metastable  population.  The 
function,  £x, 


2  PrTi 


tX=- 


(28) 


is  the  remaining  unknown  quantity  that  depends  on  an  un¬ 
known  excitation  rate  coefficient,  kxem.  The  metastable  line 
excitation  rate  coefficient  can  be  estimated  from 


k)m=PJi<  (29) 

where  PK  is  the  probability  that  the  2 p,  state  decays  through 
radiation  of  wavelength  X.  If  the  radiation  is  associated  with 
a  transition  to  the  Ly5  metastable  state,  Pk=Pmi.  Alterna¬ 
tively,  for  transitions  from  a  2 p,  to  the  \sA  state,  Px=Pri. 
Because  the  2 p{  levels  coupled  to  the  1^5  metastable  state 
have  similar  energies,  and  the  energy  intervals  between  them 
and  the  1  s5  level  are  only  approximately  1  eV,  and  because 
all  excitation  transitions  are  dipole  allowed,  we  apply  follow¬ 
ing  approximation  for  the  excitation  rate  coefficients,  Th 


TABLE  I.  Parameters  for  2/?,  states  optically  coupled  to  the  1  s5  metastable 
level  (see  Fig.  1).  In  case  of  the  two  transitions  associated  with  the  2 p7  level, 
the  840.9  and  916.3  nm  lines  arc  associated  with  the  1$5  and  ls4  lower 
states,  respectively.  Px  is  the  branching  ratio  for  radiative  decay  of  the 
respective  2 p,  level  at  wavelength,  X. 


2p6 

2p1 

2  Pfi 

2p9 

2/>io 

2^+1 

5 

3 

7 

5 

3 

0.753 

0.099  0.901 

1.0 

0.379 

0.976 

X  (nm) 

823.2 

840.9  916.3 

881.9 

904.5 

980.0 

1.890 

23.96  2.632 

1.016 

3.754 

2.430 

II 

"Ln 

2  PrTi 

i 

2  Pr,(2J,+  0 

i 

(31) 

PJj 

Pk(2Jj+\) 

The  branching  probabilities,  Px,  for  the  NIR  lines  associated 
with  2 Pi  levels  that  are  coupled  to  the  1j5  level  are  shown  in 
Table  I.  They  were  derived  from  the  emission  cross  sections 
reported  by  Fons  and  Lin.29  The  resulting  £x  values  are  also 
listed. 


C.  Sample  NIR  model  calculations 

Equation  (27)  permits  a  simulation  of  the  Xel  NIR  spec¬ 
trum  as  a  function  of  Te  and  a.  The  cross  sections  reported 
for  energies  ranging  from  10  to  70  eV  in  Part  I  are  used. 
In  order  to  extend  the  model  to  higher  electron  temperatures 
(>20  eV),  we  have  taken  the  values  reported  by  Fons  and 
Lin  for  electron  energies,  £,  between  70  and  150  eV  at  a 
nominal  pressure  p=  1  mTorr  and  modeled  them  with  the 
simple  exponential  (units  in  eV,  mTorr), 

o-x(£>70,/>=  1.0)  =  Ox(70,1.0)(^)  \  (32) 


where  the  power  cx  is  given  in  Table  II.  The  zero-pressure 
cross  sections  are  then  obtained  from  the  semiempirical  ex¬ 
pression  derived  in  Part  I, 


crk(E  >  70,0)  = 


_ crx(70,1.0)(E/70)c* _ 

{ 1  +  <jx[(E  -  1 1 .5 IE  -  9.8)]"*(  1  -  e'*)2} ' 


(33) 


ax,  /ix,  and  yx  are  line-specific  parameters  provided  in  Part  I. 
This  can  be  regarded  as  a  crude  estimate  and  a  detailed 
analysis  of  the  pressure  dependence  at  the  respective  high 
electron  energies  is  warranted  to  improve  the  model  reliabil¬ 
ity  at  high  values  of  Te. 

The  Kk(Te,a)/ fx  term  is  plotted  in  Fig.  2  as  a  function 
of  Tg  for  the  selected  NIR  lines  that  are  associated  with  the 
upper  states  that  are  optically  coupled  to  the  1^5  state.  The 
fractional  increases  are  calculated  for  a  typical  value  of  a 
=0.8.  Because  the  Xe2+  contribution  to  NIR  emissions  is 


7i  «  2Jj  +  1 ,  (30) 

signifying  that  the  transition  rate  coefficients  are  proportional 
to  the  degeneracy  of  the  upper  2 p,  level  with  total  angular 
momentum  quantum  number,  Jh  Based  on  this  assumption, 
the  function  is  given  by 


TABLE  II.  Exponential  decay  power,  cK%  for  the  selected  NIR  lines. 


X 

(nm) 

788.7 

823.2 

828.0 

834.7  881.9 

904.5 

916.3 

980.0 

-0.607 

-0.77 

-0.589 

-0.892  -1.07 

-0.7 

-0.473 

-0.86 
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FIG.  2.  The  electron  temperature  dependence  of  the  ATX/ £x  factor  for  the  five 
lines  that  are  optically  coupled  with  the  1j5  level. 

much  weaker  than  the  Xe+  contribution,  and  because  a  ap¬ 
pears  both  in  the  numerator  and  the  denominator  in  similar 
form,  Kx  is  only  weakly  dependent  on  a  [see  Eq.  (21)].  From 
the  graph  it  can  be  deduced  that  the  fractional  increase  in 
emission  intensity  due  to  the  inclusion  of  stepwise  excitation 
via  metastables  ranges  between  — ■  1 50%  and  300%. 

The  calculated  relative  line  intensities  are  also  weakly 
dependent  on  o',  as  is  demonstrated  in  Fig.  3,  where  we  com¬ 
pare  for  a=0.5  and  0.8  the  emission  rate  ratios,  7X \U\2'  be- 


FIG.  3.  Calculated  electron  temperature  dependence  of  the  intensity  ratios 
of  the  Xel  834.7,  828.0,  and  823.2  nm  lines  for  a=0.8  and  0.5.  The  hori¬ 
zontal  lines  at  low  Te  in  the  bottom  plot  are  the  pure  ionic  emission  ratios 
for  the  respective  charge  ratio  parameters.  The  ratios  of  the  bottom  plot 
based  on  the  present  stepwise  model  are  not  applicable  at  low  temperatures 
(<2  eV). 


tween  the  834.7  and  828.0  nm  lines,  which  are  uncoupled 
with  the  \s5  metastable  state  (i.e.,  Kx  =  0),  and  the  823.2  and 
828.0  nm  lines,  where  Tf823  #  0.  There  is  a  slight  a  depen¬ 
dence  up  to  5  eV,  which  is  the  temperature  region  where  ion 
collisions  produce  significant  emission  rates.  The  ratios  arc 
strongly  dependent  on  Te  over  a  significant  temperature 
range,  providing  an  accurate  diagnostic  between  1  and 
10  eV.  Below  1  eV,  the  ion  contributions  dominate,  thus  ex¬ 
plaining  the  nearly  temperature  independent  emission  ratio 
between  the  834.7  and  828.0  nm  lines.  The  823.2/828.0  line 
ratio  plot  is  presented  from  2  eV.  At  lower  temperatures,  Eq. 
(25)  is  no  longer  accurate  and  the  823.2  line  cannot  be  mod¬ 
eled  with  Eq.  (27).  At  very  low  temperatures  below  1  eV, 
where  ion  collisions  dominate  over  other  excitation  mecha¬ 
nisms,  this  ratio  approaches  the  constant  value  of  7.6  and  7.9 
for  o-=0.8  and  o,=0.5,  respectively.  These  values  are  indi¬ 
cated  by  horizontal  lines  in  Fig.  3  at  low  Te.  The  high 
823/828  ratio  at  low  Te  derives  from  the  significant  Xe+ 
emission  cross  sections  at  300  eV  for  the  823.2  nm  line  ra¬ 
diance  in  comparison  to  the  828.0  nm  line.  Above  8  eV,  the 
ion-impact  contribution  is  negligible  given  the  much  larger 
fraction  of  electrons  with  energies  above  the  excitation 
thresholds.  At  these  energies,  the  ratios  vary  only  slowly 
with  Tc  which  can  be  attributed  to  the  weaker  energy  depen¬ 
dence  of  electron  emission  cross  sections  at  higher  electron 
energies. 

A  computer  program  was  written  that  uses  a  least- 
squares  approach  to  adjust  Te ,  a  pressure  parameter  signify¬ 
ing  optical  opacity  (see  Part  I),  and  a  to  optimally  reproduce 
a  normalized  set  of  the  observed  selected  eight  NIR  line 
intensities.  The  performance  of  this  model  will  be  demon¬ 
strated  on  the  D-55  anode  layer  thruster  (TAL)  data  in 
Sec.  IV. 


D.  Visible  line  diagnostics 

The  VIS  plume  spectra  reported  in  Part  I  show  promi¬ 
nent  lines  at  467. 1  and  484.4  nm  that  are  associated  with  the 
highest  Xe+  and  Xe++  emission  excitation  cross  sections. 
These  lines  can,  therefore,  be  considered  to  diagnose  a. 
Simple  kinetic  estimates,  however,  lead  to  the  conclusion 
that  Xell  line  excitation  in  the  VIS  spectrum  must  have  very 
important  contributions  from  process  (3),  even  at  relatively 
low  electron  temperatures.  Thus,  a  proper  analysis  would 
require  knowledge  of  the  respective  electron-ion  emission 
excitation  cross  sections.  Furthermore,  the  484.4  nm  line  lies 
very  close  to  a  6s-lp  Xel  line  at  484.3  nm.  Consequently,  a 
high-resolution  experiment  is  necessary  for  an  accurate 
analysis. 

Despite  the  fact  that  cross  sections  for  electron-ion  exci¬ 
tation  processes  are  not  accurately  known,  we  can  state  that 
the  emission  excitation  rates  that  depend  on  the  electron  tem¬ 
perature  should  remain  constant  in  regions  of  electron  gas 
isothermicity.  Therefore,  Xell  line  emission  intensities 
should  be  proportional  to  the  number  densities  of  the  collid¬ 
ing  species; 
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A-N^-aN2'.  (34) 

Hence,  disregarding  the  variation  of  a  in  the  regions,  where 
the  electron  temperature  is  almost  constant,  the  charged  par¬ 
ticle  number  density  distribution  may  be  derived  from  the 
Xell  emission  intensity  measurement  as,  Ne~N\  ~ 

III.  EXPERIMENT 

Optical  spectra  of  the  TAL  were  recorded  at  the  TSNII- 
MASH  test  facility.  The  test  facility  has  a  1.8  m  diameter 
and  is  4  m  long.  It  was  pumped  by  five  oil  diffusion  pumps, 
which  resulted  in  a  base  pressure  of  1-2  X  10-6  Torr.  During 
the  operation  of  the  TAL,  the  base  pressure  grew  to  2 
X  10"5  Torr  (Xe  reference).  A  more  detailed  description  of 
the  facility  has  been  reported  previously.33  The  vacuum 
chamber  was  equipped  with  a  big  window  for  optimal  data 
acquisition  in  the  TAL  near-field  plume  and  discharge  re¬ 
gions.  General  features  of  the  TALs  as  well  as  detailed  in¬ 
formation  about  their  design  and  performance  characteristics 
may  be  found  in  the  literature.1 1_  4 

The  spectra  discussed  in  this  particular  work  were  re¬ 
corded  with  two  instruments.  The  first  was  the  asymmetric 
Ebert-Fastie  0.5  m  F/6  scanning  monochromator  equipped 
with  a  commercial  photomultiplier  tube  that  was  described  in 
Part  I.  This  instrument  was  primarily  used  for  qualitative 
analysis  of  the  TAL  spectra  and  comparison  of  the  spectra 
with  the  beam  experiment.  The  entrance  and  exit  slits  of  the 
monochromator  were  set  between  0.05  and  0.4  mm  to  opti¬ 
mally  meet  the  requirements  of  high  signal-to-noise  ratio  and 
satisfactory  spectral  resolution.  Thus,  the  slit  width  was  wid¬ 
ened  for  the  weaker  plume  region  compared  with  the  high 
intensity  discharge  region.  The  spectral  resolution  of  the 
monochromator  varied  from  about  0. 1  nm  in  the  discharge  to 
1  nm  in  the  plume  region. 

The  quantitative  studies  presented  here  were  earned  out 
with  a  charge-coupled  device  (CCD)-based  spectrometer.  A 
S2000  Ocean  Optics  spectrometer  with  600  groove /mm 
grating  was  used.  The  use  of  a  25  mm  slit  results  in  a  spec¬ 
tral  resolution  of  about  1.2  nm.  The  wide  dynamic  range  of 
the  CCD  (~105)  allowed  utilization  of  the  spectrometer  both 
for  discharge  and  plume  measurements.  The  entire  radiation 
spectra  were  measured  instantaneously  along  spans  (chords) 
of  a  selected  cross  section  perpendicular  to  the  TAL  plume 
axis.  The  spectrograph  of  the  analyzer  acquired  radiation 
through  a  fiber  optic  cable  equipped  with  a  lens  at  its  sensing 
end.  The  lens  could  be  translated  in  the  vertical  direction  by 
a  precise  stepper  motor  mechanism  while  observing  the  ra¬ 
diation  along  horizontal  lines  of  sight,  as  illustrated  in  Fig.  4. 
The  optical  system  ensured  about  a  2  mm  spatial  resolution 
in  the  discharge  and  near-field  regions  and  5-7  mm  in  the 
plume  regions,  respectively.  The  plume  region  was  located 
100-200  mm  downstream  from  the  thruster  exit  plane.  Inte¬ 
gral  (span  ward)  radiation  of  selected  spectral  lines  of  inter¬ 
est  along  different  lines  of  sight  could  be  easily  extracted 
from  these  measurements.  Radial  dependences  of  the  specific 
radiation  for  selected  lines  were  reconstructed  from  the  span- 
ward  profiles  with  the  use  of  the  standard  Abel  deconvolu¬ 
tion  procedure.4  The  wavelength  and  spectral  sensitivity  of 


Scan  Direction 

FIG.  4.  Experimental  measurement  geometry.  The  scan  direction  indicates 
the  motion  of  the  fiber  optical  sensor  cable  when  radial  profiles  are 
measured. 

the  optical  multichannel  analyzer  were  determined  in  the 
same  way  as  described  for  the  monochromator  instrument  in 
Part  I. 

IV.  RESULTS 

A.  Electron  temperature  measurements 

As  mentioned  in  Sec.  II,  the  most  appropriate  diagnos¬ 
tics  would  use  emission  lines  uncoupled  from  metastable 
states  to  avoid  errors  introduced  by  approximations  of  the 
metastable  state  excitation  functions.  In  the  case  of  the  NIR 
lines,  however,  the  uncoupled  emission  lines  at  788.7  and 
834.7  nm  exhibit  relatively  low  excitation  cross  sections, 
making  measurements  problematic  in  the  plume  region 
where  plasma  radiation  is  weak.  In  order  to  test  the  model  for 
consistency,  we  compare  electron  temperatures  derived  using 
Eq.  (18)  from  the  observed  ratio,  A34/A28’  anc*  using  Eq.  (27) 
from  the  observed  ratio,  ^23/ ^828-  Both  intensity  ratios 
/ 834/ ^828  and  were  measured  in  the  TAL  D-55  dis¬ 

charge  and  along  the  plume  axis.  The  thruster  was  operated 
at  an  anode  voltage  of  300  V,  and  a  discharge  current  of 
3.2  A.  The  electron  temperature  was  derived  from  the  mod¬ 
eled  intensity  ratios  shown  in  Fig.  3  for  a=0.8  and  assuming 
optically  thin  conditions.  Results  of  these  measurements  are 
tabulated  in  Table  III  and  plotted  in  Fig.  5.  Electron  tempera¬ 
ture  measurements  from  a  Langmuir  probe32  are  also  shown 
at  larger  axial  distances.  The  window  position  in  the  TAL  test 
chamber  precluded  measurements  at  axial  distances  exceed¬ 
ing  200  mm,  thus  preventing  a  direct  comparison  of  the 
probe  measurements,  which  were  taken  at  300  mm  and 
higher. 

All  of  the  temperatures  derived  from  the  ratio 

are  lower  than  the  corresponding  values  from  the 
ratio.  The  sources  of  the  discrepancy  will  be  discussed  in 
Sec.  V.  At  5  mm,  where  the  discharge  of  the  TAL  is  probed 
at  a  radial  distance  corresponding  to  the  annulus,  a  high  elec¬ 
tron  temperature  exceeding  50  eV  is  derived,  while  at  larger 
axial  distances,  the  electron  temperature  drops  to  values  be¬ 
tween  2.9  and  5  eV  along  the  axis  of  the  thruster.  Except  in 
the  discharge,  where  the  electron  temperature  is  high,  and, 
therefore  the  two-line  approach  is  less  accurate,  the  agree¬ 
ment  between  the  two  derivations  is  good,  validating  the 
present  treatment  of  effects  due  to  metastable  stepwise  exci¬ 
tation  mechanisms.  In  the  far  field,  the  agreement  with  the 
Langmuir  probe  measurements  is  also  very  satisfactory. 

The  pressure  dependence  of  the  two-line  approach  is 
tested  by  comparing  the  derived  temperatures  assuming  a 


Downloaded  27  Jun  2006  to  146.153.144.35.  Redistribution  subject  to  AIP  license  or  copyright,  see  http://jap.aip.org/jap/copyright.jsp 


113305-8  Karabadzhak,  Chiu,  and  Dressier 


J.  Appl.  Phys.  99,  113305  (2006) 


TABLE  III.  Measurements  of  intensity  ratios  and  derived  electron  temperature  for  a  D-55  TAL,  operated  at  an 
anode  voltage  of  300  V,  and  a  current  of  /= 3.2  A.  The  25  mm  radial  position  corresponds  to  the  thruster 
annulus. 


Position 

TAL 

T,  (eV)/CRM 

Axial 

Radial 

(mm) 

(mm) 

W  ^828 

^828 

/«34/ ^828 

his/ 1*2* 

eight  line 

5 

25 

0.231 

1.125 

50 

>70 

>70 

25 

25 

0.311 

3.35 

4.86 

5.74 

6.1  ±0.9 

25 

0 

0.332 

3.89 

3.63 

4.58 

4.9±0.7 

50 

0 

0.346 

4.48 

3.13 

3.83 

4.29±0.53 

100 

0 

0.354 

4.64 

2.94 

3.68 

3.86±0.42 

200 

0 

0.355 

4.74 

2.92 

3.61 

3.39±0.32 

nominal  pressure  of  1.5  mTorr  corresponding  to  the  optical 
conditions  of  the  beam  experiment  in  Part  I.  Not  surprisingly, 
there  is  almost  no  effect  in  the  plume  where  the  electron 
temperatures  are  insufficient  to  cause  excitation  at  energies 
where  cascading  and  the  resulting  radiation  trapping  effects 
are  significant.  Furthermore,  ion  collisions  are  very  impor¬ 
tant  at  low  temperatures,  and  they  do  not  exhibit  radiation¬ 
trapping  effects  (see  Part  I).  In  the  discharge  at  5  mm,  how¬ 
ever,  the  derived  electron  temperature  from  the  834/828 
ratio  is  lowered  to  13  eV,  while  that  for  the  823/828  ratio 
descends  to  50  eV  at  1.5  mTorr.  Since  the  dense  discharge 
region  is  most  likely  to  exhibit  signs  of  optical  opacity,  ra¬ 
diation  trapping  effects  must  be  considered  an  additional 
source  of  error  for  the  discharge. 

Given  the  errors  in  intensity  measurement  and  in  the 
measured  emission  cross  sections  of  Part  I,  and  the  approxi¬ 
mations  associated  with  the  present  pressure  dependence 
(Part  I)  and  metastable  contribution  model,  the  accuracy  of 
Te  is  expected  to  increase  with  the  number  of  lines  included 
in  the  analysis.  Figure  6  shows  the  results  of  an  eight-line 
analysis  assuming  optically  thin  conditions  on  selected  NIR 
lines  associated  with  different  2 pv  levels.  The  measurements 
are  taken  at  axial  distances  of  5  and  200  mm  (discharge  and 
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FIG.  5.  Electron  temperatures  along  the  D-55  thruster  axis  reduced  from 
measurements  of  the  respective  Xel  line  intensity  ratios,  at  an  operating 
condition  of  V=300  V,  7=3.2  A.  The  optical  diagnostic  data,  assuming 
«=0.8,  are  compared  with  the  Langmuir  probe  measurements  at  large  axial 
distances. 


plume,  respectively).  The  measured  lines,  /™L,  and  those 
corresponding  to  the  best-fit  electron  temperature,  /fM,  are 
normalized  to  l.  In  both  cases,  an  ion  charge  parameter,  a 
=0.8,  and  optically  thin  conditions  are  assumed. 

The  calculation  of  a  matrix,  (7^,a,*2),  where  \2  is  the 
relative  average  least  squares  given  by 


2  i "  (/7AL  -  /f RM)2 
x  »w  (fJAL)2  ’ 


(35) 


and  n  is  the  number  of  lines  in  the  analysis,  shows  that  Tc 
is  only  weakly  dependent  on  a.  A  temperature  of 
3.39±0.32  eV  minimizes  \2  in  the  plume  (200  mm  axial  dis¬ 
tance).  The  errors  are  the  statistical  90%  confidence  limits 
assuming  a  standard  deviation  of  10%  for  all  of  the  measured 
line  intensities,  a  highly  conservative  intensity  accuracy. 
Variation  of  a  results  only  in  minor  differences  (<6%) 
when  varied  from  0.5  to  1.0  for  the  plume  measurement.  In 
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FIG.  6.  Comparison  between  observed  and  modeled  NIR  line  intensities  for 
discharge  (5  mm  axial  distance)  and  plume  (200  mm  axial  distance)  radi¬ 
ances.  The  intensities  of  the  measured  and  modeled  spectra  are  normalized 
to  1.  The  associated  electron  temperatures,  determined  assuming  a=0.8,  are 
the  result  of  a  least-squares  fit  for  the  plume  and  are  fixed  at  40  eV,  for  the 
discharge,  where  no  least-squares  minimum  is  found  below  70  eV. 
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FIG.  7.  Standard  deviation  and  ion  fractional  radiance  contribution  as  a 
function  of  electron  temperature  determined  from  the  eight-line  NIR  model 
calculations  of  the  plume  and  discharge  line  intensities. 

the  discharge  (5  mm  distance  from  thruster  plane  outside  an¬ 
nulus),  however,  no  \2  minimum  is  found  between  1  and 
70  eV.  We  compare  the  discharge  measurement  in  Fig.  6 
with  a  set  of  spectral  intensities  calculated  for  7^=40  eV. 
At  this  energy,  99.8%  of  the  Maxwellian  electron  energy 
distribution  is  below  300  eV,  the  maximum  energy  for  which 
we  have  incorporated  cross  sections  in  the  model. 

The  sensitivity  of  the  eight-line  approach  is  demon¬ 
strated  in  Fig.  7,  where  the  relative  standard  deviation,  \ ,  is 
plotted  versus  Te  for  the  two  measurements.  The  figure  also 
plots  the  fraction  of  the  total  radiance  attributable  to  ion 
collisions.  While  there  is  not  a  \  minimum  in  the  applicable 
electron  temperature  range  for  the  discharge,  it  is  well  de¬ 
fined  for  the  plume.  At  the  derived  plume  Te  of  3.39  eV,  22% 

TABLE  IV.  Te  (eV)  determinations  in  the  TAL  plume  (200  mm  axial  dis¬ 
tance)  and  discharge  (5  mm  axial  distance)  using  the  present  CRM  model, 
“eight  lines”  corresponds  to  a  fit  to  the  experimental  eight  NIR  lines  listed  in 
Table  II  using  the  full  CRM  model  including  ion  and  metastable  contribu¬ 
tions.  “Three  lines”  refers  to  the  full  CRM  model  that  only  includes  the 
emission  transitions  that  are  not  coupled  to  the  metastable  1  s5  state.  “No 
ions”  is  a  model  that  does  not  include  ion  collision  contributions.  “No  meta” 
identifies  the  results  if  effects  due  to  metastables  are  neglected  (tfx=0).  The 
errors  are  the  statistical  90%  confidence  limits  assuming  a  10%  standard 
deviation  in  the  TAL  data  points. 


Te 

Tt 

T, 

Te 

(eight  lines) 

(three  lines) 

(no  ions) 

(no  meta) 

Plume 

3.39±0.32 

2.74±0.06 

2.49±0.40 

1 .61  ±0.21 

Discharge 

>70 

40±30 

>70 

10.5±2.2 

of  the  radiance  from  the  selected  eight-line  analysis  can  be 
attributed  to  ion  collisions.  While  the  ion  contribution  drops 
as  expected  with  Te  near  the  plume  temperature,  at  higher 
electron  temperatures  the  calculated  ion  contribution  is  seen 
in  the  top  panel  to  increase  again  after  reaching  a  minimum 
of  3%  at  20  eV.  This  can  be  attributed  to  the  decreasing 
emission  cross  sections  with  energy  as  Te  increases,  and  the 
decreasing  fraction  of  electrons  in  an  energy  range  where  the 
emission  excitation  cross  sections  are  very  high. 

The  results  of  the  TAL  evaluations  are  summarized  in 
Table  IV  where  calculations  for  an  eight-line  analysis  are 
compared  with  a  three-line  analysis  involving  only  the  2px, 
2p3,  and  2 p5  emission  lines  that  are  not  coupled  to  the  1^5 
metastable  state,  and  eight-line  models  that  neglect  ion  con¬ 
tributions  (no  ions)  and  metastables  (no  meta).  In  the  case  of 
the  plume  measurement,  the  three-line  model  produces  a 
slightly  lower  value  of  Te  than  the  eight-line  analysis.  Varia¬ 
tion  of  the  effective  pressure  between  0.1  and  1.5  mTorr 
leads  to  negligible  changes  in  the  derived  electron  tempera¬ 
ture.  Again  this  can  be  attributed  to  the  fact  that  at  low  elec¬ 
tron  temperatures,  the  electron-impact  and  ion-impact  radia¬ 
tion  processes  are  not  strongly  affected  by  radiation  trapping. 
The  determined  temperatures  from  eight-line  analyses  that 
neglect  ion  collisions  and  metastables  are  significantly  lower 
for  the  plume  models.  The  effect  of  neglecting  metastables  is 
particularly  large.  Neglecting  both  metastables  and  ions  leads 
to  a  further  reduction  in  the  derived  electron  temperature  in 
the  eight-line  analysis. 

Similarly,  neglecting  metastable  contributions  when 
modeling  the  discharge  spectrum  results  in  the  determination 
of  a  much  lower  discharge  temperature  of  ~10  eV.  Neglect¬ 
ing  ions  does  not  affect  the  analysis  due  to  the  small  ion 
contribution  at  high  Te.  As  in  the  plume  case,  a  three-line 
analysis  produces  a  lower  value  of  Te,  however,  the  accuracy 
of  the  determination,  demonstrated  by  the  large  error  of 
±30  eV,  is  very  low.  Assuming  a  higher  effective  pressure  of 
1.5  mTorr  results  in  a  lower  temperature  of  51±24eV, 
which  is  similar  to  the  observation  made  in  the  two-line 
analysis  above. 

B.  Plasma  density  measurements 

Figure  8  shows  a  radial  profile  of  the  ratio,  mea¬ 

sured  at  a  distance  of  5  mm  from  the  D-55  thruster  exit 
plane.  Interestingly,  the  ratio  is  almost  constant  across  the 
discharge  where  one  may  assume  that  the  electron  tempera¬ 
ture  stays  constant.  Isothermicity  of  the  electron  gas  in  the 
near  discharge  region  that  is  limited  by  a  20-40  mm  radial 
distance  (the  annulus  position  is  schematically  shown  by 
dashed  vertical  lines)  was  also  demonstrated  by  measure¬ 
ments  of  the  electron  temperature  profile  covering  a  10  mm 
section  by  Domonkos  et  al. 14  At  this  condition  Eq.  (34)  may 
be  applied.  Figure  9  shows  a  charged-particle  radial  density 
profile  determined  from  this  equation  at  5  mm  from  the  D-55 
TAL  (the  dashed  vertical  lines  again  point  out  the  annulus 
position).  The  Xe  ionic  line  Xell  541.9  nm  was  used  for 
these  measurements.  Utilization  of  other  Xell  lines  from  the 
VIS  spectral  region  results  in  very  similar  profiles.  The  maxi¬ 
mum  of  the  charge  density  on  the  radial  profile  is  slightly 
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FIG.  8.  Line  intensity  ratio  measurements,  lt2 across  the  D-55  TAL 
discharge  (5  mm  from  the  thruster  exit).  The  thruster  was  operated  at 
V=300  V,  7=3.2  A.  The  radial  position  of  the  annulus  is  indicated  with 
vertical  dashed  lines. 

shifted  from  the  annulus  center  towards  the  thruster  axis.  A 
similar  observation  has  been  reported  by  Haas  and  Gallimore 
in  an  SPT  using  a  fast  Langmuir  probe.34 

V.  DISCUSSION 

The  emission  excitation  cross  sections  reported  in  Part  I 
have  been  incoiporated  into  a  collisional  radiative  model 
covering  the  prominent  NIR  Xel  lines.  These  lines  are  asso¬ 
ciated  with  emission  rates  that  exceed  3  X  106  s"1  and,  there¬ 
fore,  the  model  can  neglect  depopulation  processes  other 
than  radiative  decay.  This  CRM  for  the  Hall  thruster  plasma 
radiation  diagnostics  includes  both  ion  collisional  excitation 
mechanisms  as  well  as  an  approach  to  account  for  the  effect 
of  stepwise  excitation  mechanisms  involving  metastables. 
The  model  is  applied  to  sample  spectra  observed  from  the 
discharge  and  plume  of  a  D-55  TAL  thruster  operating  in  a 
ground-based  test  facility. 

A  least  squares  fit  of  the  electron  temperature  at  fixed  ion 
charge  ratios  and  pressures  leads  to  satisfactory  reproduction 
of  the  selected  NIR  line  intensities  observed  at  various  loca¬ 
tions  outside  of  the  thruster  exit  plane.  The  derivation  at  low 
electron  temperatures  appears  only  weakly  dependent  on  the 
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FIG.  9.  Radial  distribution  of  the  charge  density  at  an  axial  distance  of 
5  mm  as  measured  from  the  541.9  nm  Xell  line.  The  D-55  TAL  was  oper¬ 
ated  at  V=300  V,  /=3.15  A.  The  radial  position  of  the  annulus  is  indicated 
with  vertical  dashed  lines. 


unknown  charge  ratio,  or,  and  radiation  trapping  effects  are 
not  significant  provided  zero-pressure  extrapolated  cross  sec¬ 
tions,  as  reported  in  Part  I,  are  used.  The  model  extracts 
electron  temperatures  that  are  consistent  with  values  typi¬ 
cally  quoted  for  such  engines,  in  particular,  for  the  near-field 
plume  region. 11,14,35  The  electron  temperatures  also  agree 
with  Langmuir  probe  measurements  determined  for  the  same 
thruster  at  the  same  operating  conditions.  Similar  values  of 
Te~ 2.6  eV  were  previously  measured  at  an  axial  position  of 
300  mm  from  the  exit  plane  of  the  D-55  thruster  using  a 
Langmuir  probe  by  Karabadzhak  et  al.} 5  Domonkos  et  al.1 4 
measured  radial  profiles  of  electron  temperature  at  several 
axial  distances  from  a  D-55  thruster.  At  20  mm,  they  report 
an  axial  electron  temperature  of  3.5  eV,  while  at  a  radial 
distance  corresponding  to  the  annulus,  the  temperature 
climbs  to  5.7  eV.  These  measurements  are  comparable  with 
present  values  of  4.9 ±0.7  and  6.0±0.9  eV,  respectively,  de¬ 
termined  in  the  eight-line  analysis  at  25  mm  (Table  III). 
Similar  temperatures  have  been  modeled  by  Karabadzhak 
et  al ,35  and  were  also  used  as  input  for  successful  computa¬ 
tions  of  the  D-55  thruster  plume  parameters  by  DSMC-PIC 
methods.17 

Closer  to  the  thruster,  probe  measurements11  14  result  in 
electron  temperatures  exceeding  10  eV.  Gamer  et  al.  11  re¬ 
ports  that  the  temperature  increases  up  to  80  eV  and  higher 
in  the  discharge,  again  consistent  with  the  present  model  and 
observations.  High  TAL  discharge  temperatures  have  also 
been  predicted  by  Choueiri.36  Domonkos  et  al.] 4  reported  a 
non-Max wellian  behavior  of  the  probe  current  in  this  near¬ 
field  region.  At  such  conditions  it  can  be  expected  that  probe 
and  optical  measurements  produce  different  results  since 
probes  tend  to  be  more  sensitive  to  the  energy  distribution  of 
low-energy  electrons,  while  the  optical  measurements  pro¬ 
vide  a  measure  of  the  electron  energy  distribution  above 
threshold.  From  inspection  of  the  two-line  intensity  ratios 
shown  in  Fig.  3,  it  can  be  concluded  that  the  intensity  ratios 
predicted  by  the  NIR  model  become  insensitive  to  electron 
temperature  variations  beyond  15-20  eV.  A  more  accurate 
model  of  higher  electron  temperatures  would  depend  on  an 
analysis  of  visible  emission  lines,  including  Xell  lines.  A 
model  including  the  latter  depends  on  emission  cross  sec¬ 
tions  for  electron  impact  excitation  of  Xe\  which  are  not 
available. 

As  seen  in  Fig.  2,  inclusion  of  metastables  results  in 
~300%  increase  in  the  modeled  emission  intensity  of  the 
823.2  nm  line.  Comparison  of  the  present  model  with  calcu¬ 
lations  that  either  neglect  ion  contributions  or  stepwise  meta¬ 
stable  excitation  exhibit  marked  differences  in  the  derived 
electron  temperature,  generally  leading  to  significantly  lower 
values  (see  Table  IV),  if  not  complete  failure  to  reproduce 
the  spectrum.  This  work,  in  conjunction  with  our  earlier 
report,9  establishes  the  importance  of  ion-atom  collisions  in  a 
Hall  thruster  radiative  model.  The  present  model  suggests 
that  —30%  of  the  NIR  radiance  is  due  to  ions  at  a  typical 
plume  temperature  of  3  eV.  The  ion-collision  contributions 
are  a  direct  consequence  of  the  large  Xe++Xe  NIR  emission 
excitation  cross  sections  reported  in  Part  I,  and  the  signifi¬ 
cant  differences  observed  in  the  ion  and  electron-induced 
spectral  intensity  distributions  due  to  the  different  excitation 
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FIG.  10.  Correction  factor,  for  uncoupled  line,  \  =  834.7  nm,  that 

brings  Tt  determined  from  the  Igy^/I^s  ratio  m  agreement  with  that  deter¬ 
mined  from 

mechanisms.  Consequently,  intensity  ratios  between  particu¬ 
lar  line  pairs  change  dramatically  as  Te  is  increased  from  low 
values,  where  ion  impact  dominates  the  spectrum,  to  values 
where  electron  impact  is  the  primary  source  of  NIR  radiance. 

The  discrepancy  between  the  electron  temperatures  de¬ 
rived  using  the  anc*  ^834^828  intensity  ratios  is  symp¬ 

tomatic  of  the  fact  that  the  current  model  is  still  incomplete. 
Likely  sources  of  the  discrepancy  are  the  assumption  of  Eq. 
(30),  which  provides  an  estimate  for  the  relative  line  excita¬ 
tion  cross  sections  from  the  \s5  metastable  state,  and  the 
approximation  that  those  lines  that  are  not  coupled  to  the  1  s5 
level  do  not  experience  an  enhancement  due  to  some  step¬ 
wise  excitation  mechanism.  While  there  have  been  no  mea¬ 
surements  or  calculations  of  emission  excitation  cross  sec¬ 
tions  out  of  metastable  states  of  Xe,  there  are  reasons  to 
reevaluate  the  latter  assumption  for  the  2 p{  (788.7  nm)  and 
2 p3  (834.7  nm)  lines  that  the  model  considers  uncoupled 
from  metastable  states.  These  levels  are  associated  with  the 
7=  1/2  ion  core  limit,  i.e.,  they  are  referred  to  as  6 p'  states, 
and  are  mixed  with  3 pt  levels  associated  with  the  7=3/2  ion 
core  limit,  or  3/?,(7/?)  states,  as  evidenced  by  emissions  from 
6//  to  6s  and  from  Ip  to  65 '  levels.  Consequently,  while  the 
2p j  and  2 p3  6 p'  levels  are  decoupled  from  the  153(6j') 
metastable  state,  they  are  not  fully  decoupled  from  the 
155(6a’)  metastable  state.  Thus,  the  associated  788.7  and 

834.7  nm  lines  are  expected  to  be  enhanced  through  a  step¬ 
wise  mechanism,  while  the  2p5(6p)  828.0  nm  line  is  not. 

In  order  to  assess  the  magnitude  of  the  enhancement,  the 
Kk/£x  term  [Eq.  (27)]  of  the  \  =  834.7  nm  line  is  raised  from 
zero  to  bring  the  derived  value  of  Te  based  on  the 
ratio  in  agreement  with  that  based  on  The  thus  de¬ 

rived  /^834/£834  values  for  the  /8 34/^28  measurements  of 
Table  III  are  plotted  in  Fig.  10.  Interestingly,  a  decline  is 
observed  from  a  value  of  ~7%  at  ~3.6  eV  to  a  value  of 
3.5%  at  5.7  eV.  This  is  consistent  with  the  forbidden  nature 
of  the  ls5(6s)-2p3(6pr)  transition,  where  the  electron  en¬ 
ergy  dependence  of  the  excitation  cross  section  can  be  ex¬ 
pected  to  exhibit  a  sharp  maximum  above  threshold  with  a 
rapid  decrease  at  higher  energies  where  spin-changing  colli¬ 
sions  become  less  probable. 

The  discharge  data  point  in  Fig.  7,  representing  an  elec¬ 


tron  temperature  larger  than  70  eV  as  obtained  from  the 
7823//828  analysis,  is  not  shown.  A  mere  2%  enhancement  of 
the  /g34//828  ra^°  in  the  discharge  would  result  in  an  increase 
of  the  derived  Te  from  50  to  70  eV.  This  clearly  suggests 
that  this  ratio  is  not  suitable  for  an  analysis  at  high  electron 
temperatures.  The  minor  enhancements  required  to  bring  the 
temperatures  derived  from  the  a°d  ^823^828  ratio 

analyses  into  agreement  are  consistent  with  the  weaker  cou¬ 
pling  and  the  higher  respective  \s5-^2p-s(6p')  excitation  en¬ 
ergies,  leading  to  lower  cross  sections.  Incorporation  of  a  7% 

834.7  nm  line  enhancement  in  the  eight-line  analysis  of  the 
200  mm  plume  measurement  results  in  only  a  negligible  re¬ 
duction  in  the  derived  temperature  of  3.36±0.32  eV  com¬ 
pared  with  the  uncorrected  value  of  3.39 ±0.32  eV. 

The  enhancement  of  the  788.7  nm  line  can  be  expected 
to  be  even  less  due  to  the  lower  degeneracy  (7=0)  of  this 
state.  Arbitrarily  setting  the  enhancement  to  half  that  of  the 

834.7  nm  line  results  again  in  a  negligible  reduction  of  the 
derived  temperature.  From  the  above  it  can  be  concluded  that 
the  2pi(6pr)  lines  are  indeed  slightly  affected  by  a  stepwise 
mechanism,  but  the  correction  does  not  affect  significantly 
the  eight-line  analysis.  Overall,  given  the  weakness  of  the 

788.7  and  834.7  nm  lines,  the  error  in  the  intensity  measure¬ 
ment  can  be  regarded  comparable  to  the  enhancement  due  to 
stepwise  excitation.  The  present  TAL  measurements  have  de¬ 
tected  interference  from  an  additional  line  near  834.7  nm, 
most  likely  the  6p2D^/2—*5d2F^i2  emission  line  of  Xe  II. 
This  interference  was  not  resolved  in  the  experiments  that 
measured  the  cross  sections.  Although  the  interference  was 
taken  into  account  in  the  TAL  measurements,  it  adds  to  the 
error  of  the  derived  intensity  of  this  weak  line.  This  further 
emphasizes  the  need  for  an  analysis  based  on  several,  in¬ 
tense,  2 pj  lines  such  as  the  present  eight-line  analysis, 
thereby  mitigating  errors  in  intensity  derivations  as  well  as 
errors  in  the  cross  sections  and  associated  pressure  depen¬ 
dences  reported  in  Part  I. 

The  most  important  shortcoming  of  the  present  CRM  is 
the  inability  of  modeling  Xell  emission  intensities  in  the  VIS 
and  UV  spectra.  This  is  attributed  to  the  complete  absence  of 
accurate  calculations  or  measurements  of  e"+Xe+  emission 
excitation  cross  sections.  A  set  of  electron-impact  Xell  emis¬ 
sion  cross  sections,  similar  to  those  reported  in  Part  I,  would 
provide  a  wider  choice  of  transitions  that  should  also  im¬ 
prove  the  accuracy  at  high  electron  temperatures  due  to  the 
higher  excitation  thresholds  of  the  respective  ionic  levels. 

Nevertheless,  the  charge  density  distribution  can  be  de¬ 
rived  in  the  TAL  discharge  based  on  the  electron  gas  isother- 
micity  assumptions  of  Eq.  (34).  The  peak  of  the  distribution 
is  slightly  shifted  from  the  annulus  center  toward  the  thruster 
axis.  Such  behavior  of  the  charge  density  distribution  is  often 
observed  in  TALs.  It  demonstrates  a  tendency  to  form  a  high- 
density  plasma  region  along  the  axis  in  the  very  near  plume 
field.  This  region  is  conventionally  called  “crossover”  and  is 
readily  observed  by  the  naked  eye.  Outside  the  isothermicity 
region  the  results  may  only  have  a  qualitative  character.  As 
Fig.  8  illustrates,  the  intensity  ratio  and  electron  temperature 
is  highly  variable  closer  to  the  axis.  At  large  axial  distances 
the  approach  based  on  the  assumption  that  process  (3)  is  the 
only  contributor  is  expected  to  fail,  as  well.  Even  if  the  con- 
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dition  of  electron  gas  isothermicity  holds  at  large  distances, 
the  charge-ratio  parameter,  a,  may  vary  in  the  radial  direc¬ 
tion  due  to  the  dependence  of  the  angular  ion  intensity  dis¬ 
tributions  on  charge  state. 

The  ion-Xe  impact  may  also  produce  a  significant  frac¬ 
tion  of  Xell  radiance  at  low  electron  temperatures.  A  reduc¬ 
tion  of  ion-impact  interference  can  be  achieved  by  basing  the 
analysis  on  Xell  lines  that  have  low  ion-impact  cross  sec¬ 
tions  such  as  the  504.5  nm  line  (see  Part  I).  That  line,  how¬ 
ever,  is  barely  detected  in  the  plume  spectrum.  In  the  VIS 
region,  the  latter  spectrum  does  in  fact  closely  resemble  the 
Xe2+-impact  spectrum,  suggesting  further  that  process  (3)  is 
not  the  sole  contributor  to  Xell  lines  in  the  far  field.  As 
mentioned  in  Part  I,  a  detailed  analysis  of  these  conditions 
could  lead  to  a  CRM  for  the  visible  emissions  that  would 
result  in  the  additional  diagnostic  of  ion  charge  ratios  in  the 
plume.  This  work  is  reserved  for  future  research. 

VI.  CONCLUSION 

This  paper  represents  an  attempt  at  incorporating  a  broad 
set  of  emission  excitation  cross  sections  in  a  CRM  for  HET 
Xe  plasmas.  A  detailed  model  is  presented  for  a  set  of  intense 
NIR  emission  lines.  The  model  establishes  the  importance  of 
ion  processes  in  producing  emissions  at  plasma  electron  tem¬ 
peratures  below  5  eV.  The  model  also  takes  a  stepwise  exci¬ 
tation  via  metastables  into  account.  Due  to  the  markedly  dif¬ 
ferent  NIR  spectra  produced  by  ions  and  electrons,  the 
method  is  particularly  accurate  at  low  temperatures  where 
ionic  contributions  are  important,  and  thus  line  ratios  depend 
strongly  on  the  effective  electron  temperature.  The  diagnos¬ 
tic  presents  a  choice  of  measuring  eight  selected  NIR  lines 
associated  with  the  2/?,  (Paschen  notation)  upper  levels,  or 
just  two  lines  for  which  the  intensity  ratio  is  highly  Te  de¬ 
pendent.  An  eight-line  analysis  is  recommended  in  order  to 
overcome  inaccuracies  in  the  pressure  dependent  effective 
cross  sections  and  plasma  line  intensity  measurements.  For  a 
simple  two-line  analysis,  the  823.2  and  828.0  nm  lines  are 
recommended  due  to  their  close  spectral  proximity,  thereby 
minimizing  the  spectral  sensitivity  calibration  errors,  and  due 
to  the  significant  Te  dependence  at  low  temperatures.  The 
present  application  of  the  model  to  TAL  D-55  measurements 
yields  Te  values  that  are  within  errors  in  agreement  for  the 
recommended  two-line  and  eight-line  approaches.  It  is  also 
demonstrated  that  spatial  plasma  density  distributions  can  be 
determined  from  Xell  line  intensities  at  regions  of  high  elec¬ 
tron  temperatures. 
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